Diffusion Time-Scale of Porous Pressure-Sensitive Paint 
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Introduction 


Pressure-sensitive paint (PSP) is an optical pressure sensor that utilizes the oxygen 
quenching of luminescence. PSP measurements in unsteady aerodynamic flows require fast time 
response of the paint [1,2,3]. There are two characteristic time-scales that are related to the time 
response of PSP. One is the luminescent lifetime representing an intrinsic physical limit for the 
achievable temporal resolution of PSP. Another is the time-scale of oxygen diffusion across the 
PSP layer. When the time-scale of oxygen diffusion is much larger than the luminescent 
lifetime, the time response of PSP is controlled by oxygen diffusion. In a thin homogenous 
polymer layer where diffusion is Fickian, the oxygen concentration [0 2 ] can be described by 
the diffusion equation in one-dimension. 


d[0 2 J n d 2 fQ 2 J 

dt dz 2 


( 1 ) 


where D m is the diffusivity of oxygen mass transfer in the polymer, t is the time, and z is the 
coordinate normal to the layer. For an arbitrary time-dependent change of the oxygen 
concentration at the air-paint interface, the complete convolution-type solution of Eq. (1) with 
the suitable boundary conditions can be obtained by using the Laplace transform [4,5], The 
special trigonometrical-series-type solutions for a step change and a sinusoidal change of oxygen 
were used for PSP dynamic analysis [1,2,3], The solution of the diffusion equation gives a well- 
known square-law estimate for the diffusion time-scale T diff through a PSP layer, 

T diff x h 2 / D „, , ( 2 ) 

where h is the PSP layer thickness. The response time of a homogenous PSP is proportional to 
the square of the layer thickness and inversely proportional to the diffusivity of mass transfer D m . 
This estimate is correct for a conventional homogenous PSP. 
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Compared with a conventional PSP, a porous PSP has a much shorter diffusion response 
time ranging from 18|is to 500|is [6], The shorter diffusion time-scale is essentially related to an 
enlarged air-polymer interface in a porous PSP. Interestingly, recent measurements of the 
response time for three polymers GP197, GP197/BaS04 mixture and Poly(TMSP) showed that 
the classical square-law estimate (2) does not hold for porous PSPs [7, 8], As shown in Fig. 1, 
measurements gave the power-law relations T diff x h 183 for GP197, z diff ~h 107 for 

GP197/BaS04 mixture, and x diff °c h 029 for Poly(TMSP) at a temperature of 313. IK. For the 

GP197 silicone polymer, the power-law exponent is close to 2 as predicted by the classical 
estimate. Flowever, the power-law exponents for the porous polymer materials GP197/BaS04 
mixture and Poly(TMSP) are significantly smaller than 2. In addition. Fig. 2 shows that the 
power-law exponent for the polymer Poly(TMSP) linearly increases with temperature over a 
temperature range from 293. IK to 323. IK. Unfortunately, a comprehensive and universal theory 
of diffusion in porous materials has not been developed yet. In order to understand the time 
response of a porous PSP, nevertheless, this note attempts to derive the expressions for the 
effective diffusivity and the diffusion response time of a porous polymer layer from a standpoint 
of phenomenology. 


Effective Diffusivity 

Diffusion in a porous material can be considered as a diffusion problem in a two-phase 
system made up of one disperse phase and one continuous polymer. In PSP, the disperse phase 
is composed of numerous pores filled with air. An element of a porous polymer layer of the 
length /, width /, and thickness h is considered, as shown in Fig. 3. The coordinate z is normally 
directed to the polymer layer from the upper surface of the layer. We assume that many 
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cylindrical (tube-like) pores are distributed in the element and the pores are oriented in the z- 
direction. The effective radius and depth of a pore are denoted by r and h pore , respectively. 
The radius of a pore is much larger than the size of a molecule of oxygen. In general, the depth 
of a pore is smaller than the layer thickness, i.e., h pore < h . For simplicity of expression, the 

normal directional derivative of the oxygen concentration [0 2 ] at the air-polymer interface is 
denoted by 


v n ( z ) = 


d[0 2 ] 

d n 


(3) 


The effective diffusivity of the porous polymer layer with cylindrical pores is given by a 

balance equation between the mass transfer through the apparent homogenous upper surface and 
the total mass transfer across the air-polymer interface, 

^ meff ' 


Dmeffl Vn (0) D m (l N 7rr ) v n (0) + D m N ki pore v n (h pore ) 


+ D N 2k r 

pore pore 


C hpore 

v n (z)dz 

Jo 


(4) 


where N pore is the total number of the pores in the element and D m is the diffusivity of the 
polymer continuum. The integral term in Eq. (4) is the total mass transfer across the peripheral 
surface of the pores in the element. Thus, the effective diffusivity D mc/ j is given by 


D meff /D m = l + [ v n (h pore )/v n (0) - 1] N pore k r 2 pore l~ 2 


+ N 2n v l 

' iy pore r pore 1 v 


i f l pore 

n(0)\ V n (z)dz 

Jo 


(5) 


In a simplified case where v n (z) = constant across the thin layer, Eq. (5) becomes 


D m eff / D m = I + 2 “v r P L h. 


( 6 ) 
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where a v = N pore n r p 0re h pore l 2 h 1 is the volume fraction of the cylindrical pores in the 

polymer layer. Eq. ( 6 ) indicates an increase of the effective diffusivity that is proportional to the 
volume fraction of the pores and a ratio between the polymer layer thickness and the radius of 
the pore. The expression (6) for D mef - f is valid only for an ideal porous polymer layer with the 

straight cylindrical pores oriented normally. Nevertheless, this model should be generalized for 
real porous polymers where topology of the pores is often highly complicated. 

In a more realistic model, the topological structure of a pore is considered as a highly 
convoluted and folded tube in the polymer layer while the cross-section of the tube remains 
unchanged. The integral in Eq. (5) should be replaced by the integral along the path of the 
highly convoluted tube-like pore. In this case, the concept of the fractal dimension should be 
introduced because the length of the highly convoluted tube is no longer proportional to the 
linear length scale of the tube in the z-direction (e.g. h pore ) [9]. According to the length-area 

relation for a fractal path, the integral along the path is proportional to A p £/ e 2 or h p £ re , where 

dp. (1 < dp. < 2) is the fractal dimension of the path of the pore and A pore °c h 2 pore is the 

characteristic area covering over the path. Loosely speaking, the fractal dimension represents the 
degree of complexity of the pore pathway. In order to take the fractal nature of the pores into 
account, Eq. (5) is generalized by using a Riemann-Liouville fractional integral of the order d fr , 
i.e.,[10] 

D meff /£>,„ = l + [ v„ (h pore )/v n (0) - 1] N pore n r 2 pore r 2 

+ N pore r pore ^ V n‘ (0) £ ^ V„ ( z)(dzf fr 
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1—df 

Note that a unitary constant with the dimension [ m " ] is implicitly embedded in the third 

term in the right-hand side of Eq. (7) to make Eq. (7) dimensionally consistent. This 

dimensional constant is implicitly contained in all the results derived from Eq. (7). In a 
simplified case where v n (z) = constant across the thin layer, a generalized expression for D meJ j 
is 


Dmeff _ j + ^ a V r pore 


D. 


r(l + dj r ) h 


h d f r 1 

Pore | ^df r 


( 8 ) 


where r(l + d fr ) is the gamma function. Here, h pore is interpreted as a linear length scale of 


the convoluted tube in the z-direction and a v is the volume fraction of the apparent cylindrical 
pores. The expression (8) clearly shows that the effective diffusivity D mrjj is not only 


proportional to ti fr , but also related to the parameters of porosity a v r pore and h pore / h . For 


*pore ' 


dp. = 1 , Eq. (8) is simply reduced to Eq. (6). 


Diffusion Time-Scale 

For a porous polymer layer where diffusion is Fickian under some microscopic 
assumptions [ 1 1,12], the diffusion equation ( 1 ) is still a valid phenomenological model as long as 
the diffusivity D m is replaced by the effective diffusivity /T Jv;( . Hence, the estimate for the 
diffusion time-scale for a porous PSP layer is 


h 2 /D„ 




l + - 


2 a 


r - 1 

V ' pore 


h 


dp-i 


pore 


r(l + d fr ) h 


h d fr 


(9) 
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Eq. (9), a generalized form of Eq. (2 ), clearly illustrates how the parameters of porosity a v r p * re 
and h pore /h, and the fractal dimension d fr affect the response time of a porous PSP. For 
a v r~g re « 1 or h pore / h « 1 , Eq. (9) naturally approaches to the classical square-law estimate 
(2) for a homogenous polymer layer. On the other hand, for a v r p ° re » 1 and h pore / h~ 1 , 
another asymptotic estimate for T diJ j is a simple power-law 

~ h™' / D m . (10) 

The estimate (10) is asymptotically valid for a very porous polymer layer. The exponent in the 
power-law relation between the response time T di g and thickness h deviates from 2 by the fractal 

dimension d fr because of the presence of the fractal pores in the polymer layer. The relation 
(10) provides an explanation for the experimental finding that the exponent q in the power-law 
relation T djff °= h q is less than 2 for porous PSPs [7, 8]. Also, this relation can serve as a useful 
tool to extract the fractal dimension of the tube-like pores in a very porous polymer layer from 
measurements of the diffusion response time. The fractal dimension d fr of the pore in the 

polymer Poly(TMSP) is d g r =1.71. For GP197/BaS04 mixture, the fractal dimension d^ r is 

close to one. In addition, based on the experimental results shown in Fig. 2, one knows that the 
fractal dimension d fr for the polymer Poly(TMSP) linearly decreases with temperature in a 
temperature range from 293. IK to 323. IK. This implies that the geometric structure of the pore 
in Poly(TMSP) may be altered by temperature change. Note that the diffusivity D m of oxygen 

mass transfer is also temperature-dependent, but it is independent of the coating thickness h. 
Therefore, the experimental results in Fig. 2 mainly reflect the effect of temperature on the 
geometric structure of the pores in the polymer rather than the diffusivity. 
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Conclusions 


The Note develops a simple phenomenological model for the effective diffusivity of a 
porous PSP and gives a new expression that clearly illustrates the relationship between the 
diffusion time-scale and the fractal dimension of the pores. The theoretical results can not only 
explain why a porous PSP is able to achieve a very fast time response, but also quantitatively 
show how the fractal dimension and the parameters of porosity affect the response time. For 
very porous PSPs, the classical square-law estimate of the diffusion time-scale should be 
replaced by the generalized relation derived in the Note. 
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Captions 


Figure 1. Power-law relationship between response time and coating thickness for three 
polymers GP197, GP197/BaS04 mixture and Poly(TMSP) at a temperature of 313. IK. 
Figure 2. The exponent of the power-law relation between time-scale and coating 
thickness for the polymer Poly( TMSP ) as a function of temperature. 

Figure 3. An element of a porous polymer layer. 


10 



Time constant, msec 



Figure 1. Power-law relationship between response time and coating thickness for three 
polymers GP197, GP197/BaS04 mixture and Poly(TMSP) at a temperature of 313. IK. 
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Figure 2. The exponent of the power-law relation between time-scale and 
coating thickness for the polymer Poly(TMSP) as a function of temperature. 
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Figure 3. An element of a porous polymer layer. 
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